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We report measurements of the photoinduced change in reflectivity of an untwinned single crystal 
of Y Ba-zCuzO§.5 in the ortho II structure. The decay rate of the transient change in reflectivity is 
found to decrease rapidly with decreasing temperature and, below T c , with decreasing laser intensity. 
We interpret the decay as a process of thermalization of antinodal quasiparticles, whose rate is 
determined by an inelastic scattering rate of quasiparticlc pairs. 
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Much of current researchon cuprate superconductivity 
focuses on the low-energy electronic structure. Perhaps 
the most pressing question is whether the pseudogap is 
a precursor of the superconducting gap, or has an en- 
tirely different origin. Another key question is whether 
quasiparticles live in a two-dimensional environment, or 
one with pronounced one-dimensional properties on short 
length and time scales. 

The relaxation dynamics of non-equilibrium quasi- 
particles can be extremely sensitive to such local elec- 
tronic properties plp|. A promising technique to study 
quasiparticle relaxation is optical pump and probe spec- 
troscopy. In this type of experiment a pulse of light stim- 
ulates a non-equilibrium state with excess energy and 
quasiparticle density. A weaker pulse, which is delayed 
relative to the pump, probes the nonequilibrium state by 
detecting changes in the absorptance or reflectance. Dy- 
namics can be measured with sub-picosecond time reso- 
lution while anisotropy can be explored by varying the 
polarization of the pump and probe beams. 

In this paper we report pump and probe measure- 
ments performed on an untwinned single crystal of 
YBa,2Cu306.5 (Ortho II). In this material the charge 
reservoir layer is an ordered array of alternating filled and 
empty copper oxygen chains. Together with YBa 2 Cu40s 
it is one of two underdoped materials in which doping 
does not introduce disorder. We find, in contrast with 
an earlier report [p[ on underdoped films, that the decay 
rate of the photoinduced change in reflectivity depends 
strongly on the excitation density and temperature. Be- 
cause it vanishes in the limit that both tend to zero, we 
interpret the decay rate as an inelastic scattering rate of 
quasiparticle pairs. We present evidence that the decay 
corresponds to a process of thermalization, rather than 
recombination, of the nonequilibrium quasiparticle popu- 
lation. The dynamics of thermalization change suddenly 
upon crossing from the pseudogap regime to the super- 
conducting state. 



We measured the decay of the change in reflectivity at 
a photon energy of 1.5 eV due to excitation with pho- 
tons of the same energy. The optical pulses, produced by 
a mode-locked ThSapphire laser, have duration of 100 
fs, repetition rate 12 ns, and center wavelength 800 nm. 
Both the pump and probe beams were focused onto the 
sample with a 20 cm focal length lens, yielding a spot 
size of 75 /i diameter. A photoelastic modulator varied 
the pump intensity at 100 kHz and a vibrating mirror 
oscillating at 40 Hz varied the time delay between pump 
and probe. This double modulation provides sensitivity 
to the fractional change in reflectivity AR/R of 10~ 7 . 
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FIG. 1. (a) AR/R vs. t at constant pump intensity for 
different temperatures (polarization perpendicular to chains) . 
(b) Maximum amplitude of AR/R vs. T for both polarization 
parallel and perpendicular to the chains. 

Fig. la shows AR/R due to absorption of a pump 
pulse of energy 0.07 nJ vs. time (t), at several tem- 
peratures. These reflectivity transients were measured 
with pump and probe polarization each aligned perpen- 
dicular to the chain axis but the decay rate is the same 
when measured with polarization along the chains. The 



amplitude of AR/R increases and the decay rate de- 
creases as the temperature is lowered. Fig. lb shows 
the temperature dependence of the maximum amplitude 
(AR/R) max . The reflectivity change is smaller when 
measured with pump and probe parallel to the chains, 
particularly above T c . 

The T dependence shown in Fig. 1 is consistent with 
earlier reports of AR/R in underdoped (twinned) thin 
films excited at 1.5 eV and probed at either near- infrared 
M or far-infrared H frequencies. These data have re- 
vealed a strong correlation of AR/R with the strength of 
spin and current fluctuations at the characteristic energy 
scale for superconductivity in the cuprates. As empha- 
sized in Ref. [Q, the increase of AR/R with decreasing 
T tracks the growth of a spin resonance peak seen by 
neutron scattering [p|,p| and a gap-like feature seen in op- 
tical conductivity [|7|j8|]. Furthermore, AR/R shows the 
same behavior with varying carrier concentration as do 
these spectroscopic features, namely abrupt onset at T c 
for near optimal doping |9| and a slow onset above T c 
for underdoped materials |ljJ4| . Such considerations sug- 
gest that 1.5 eV photons ultimately generate low-energy 
excitations, despite the fact that they initially create an 
electron-hole pair whose energy is much larger than the 
maximum gap energy Ao- 

It is generally agreed that the generation of low en- 
ergy quasiparticles occurs via a rapid cascade process, in 
which the parent electron-hole pair creates a large num- 
ber of excitations at much lower energy. What has re- 
mained less clear is the mechanism by which this nonequi- 
librium population of quasiparticles causes a change in 
reflectivity in the visible region of the spectrum. How- 
ever, recent measurements B of the photoinduced change 
in far-IR reflectivity point to a mechanism for changes 
in reflectivity at 1.5 eV. Because it is crucial to under- 
standing the decay of AR/R with time, we describe this 
mechanism below. 

The far-IR measurements reveal that photoexcitation 
shifts conductivity spectral weight from low energies 
(most likely the condensate contribution at zero energy) 
to ^100 meV. This change corresponds to filling in the 
gap feature in the optical conductivity. The Kramers- 
Kronig relations dictate that this redistribution of spec- 
tral weight changes the real part of the dielectric function 
at much higher frequencies (o>) according to: 



AeiH = - 
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where Acti is the change in the real part of the optical 
conductivity. Eq. 1 shows that Aei is proportional to 
the second moment of the change in the conductivity, 
or E*- 2 ). Because the spectral weight shifts from zero to 
^100 meV, we can relate the second moment to the spec- 
tral weight shift, i.e., £ (2) w (lOOmeV) 2 J Q+ duoAa^u). 
Finally, using literature |1(| values t\ — 1.1 and 62 = 2.7 



at 1.5 eV, we can establish an approximate calibration 
between AR/R and J Q+ dujAai(uj). 

The data shown in Fig. 2 provide a test of this 
mechanism for reflectivity change at 1.5 eV. We plot 
(AR/R) max at 15 K vs. excitation density expressed 
in units of photons per planar Cu. The right-hand scale 
shows the spectral weight shift that, according to the 
preceding analysis, generates the change in reflectivity 
shown on the left-hand scale. The spectral weight is nor- 
malized to the condensate spectral weight assuming a 
penetration depth of 1600 A Q. Notice that AR/R de- 
parts from linearity at ~ 1.5 x 10 -4 , which corresponds 
to a shift of 30% of the condensate spectral weight. This 
reasonable value for the onset of nonlinearity lends cred- 
ibility to the calibration. It is important to note that 
a measurable reflectivity change occurs only if spectral 
weight shifts to frequencies of the order of Aq. If in- 
stead the spectral weight were shifted to the quasiparti- 
cle Drude peak, AR/R would be smaller by the factor 
(1/rAo) 2 , which is ~ 10 6 for these exceptionally clean 
ortho II crystals [psf. 




-a 

N 



o 



'a) 



o 

a. 



1 



Photon density (10" /planar Cu) 

FIG. 2. Ai?/_R(black circles) vs. excitation density. 
Right-hand scale shows the spectral weight shift necessary 
to produce the reflectivity change on the left. 

Given this interpretation of the change in reflectivity, 
we now turn to the central result of this paper, namely 
the rate of decay of AR/R measured as a function of T 
and pump intensity /. We have seen in Fig. 1 that the 
decay of AR/R becomes slower as T decreases. Fig. 3a, 
a plot of the reflectivity change at 15 K created by pulses 
of varying intensity, shows that the decay rate also slows 
as the excitation density decreases. 

The decays shown in Fig. 3a are non-exponential 
in that the instantaneous rate of decay, -f(t) = 
—dln(AR)/dt, decreases as t increases. To characterize 
the intensity dependence of the decay rate in a relatively 
simple way, we focus on the initial decay rate, or 7(0). 
We estimate 7(0) from the slope of a linear fit to AR/R 
between 100% and 75% of (AR/R) max , normalized by 
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FIG. 3. (a) Decay of AR/R for several pump intensities at 
15K (b)Initial decay rate 7(0) vs. (AR/R) max . 



the average value of AR(t)/R in that interval. Given our 
time resolution, this procedure provides a good estimate 
of 7(0) for rates slower than ~ lps -1 . Fig. 3b shows 7(0) 
as a function of (AR/R) max , obtained by analysis of the 
curves similar to the ones in Fig. 3a. 7(0) decreases 
linearly with (AR/ R) max and extrapolates to a nonzero 
intercept as the excitation density approaches zero. 

The dependence of 7(0) on excitation density is consis- 
tent with two-particle kinetics, for example a process in 
which a pair of quasiparticles scatter off of each other. 
For high excitation densities the rate of such scatter- 
ing is proportional to the nonequilibrium density. As 
the excitation density decreases, becoming smaller than 
the thermal population, 7(0) approaches the equilibrium 
scattering rate at that temperature |0|. 

Fig. 4 presents the intensity dependence of the de- 
cay rate for a range of temperatures, in the form of a 
logarithmic plot of 7(0) vs. T. Points with the same 
(open) symbol show 7(0) measured at the same excita- 



tion density, but at different temperatures. The data sets 
coincide at high temperatures but diverge suddenly at T c , 
showing that the dependence of 7(0) on I starts abruptly 
upon crossing from the pseudogap to the superconduct- 
ing state. The solid symbols show the decay rate in the 
limit of zero excitation density, 7^, as determined from 
the intercept of plots such as in Fig. 3b. The same quan- 
tity is plotted on linear axes in the inset. With decreasing 
temperature j t h begins to decrease in the normal state 
and is proportional to T 3 for temperatures less than ~40 
K. 
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FIG. 4. Initial decay rate 7(0) vs. T at fixed I (open sym- 
bols). Decay rate in low I limit, fth, vs. T (solid symbols in 
main panel and inset). 

Several aspects of the decay rate of AR/R suggest a 
strong connection to quasiparticle dynamics, as well as 
structure in the low-energy density of states. For exam- 
ple, the decay rate appears sensitive to the opening of 
the pseudogap, decreasing near T* in a manner which 
is remarkably similar to the nuclear spin relaxation rate 
in underdoped YBCO 14 1. The onset of superconduc- 
tivity also affects the decay dynamics — the decay rate 
depends on excitation density only below T c . Finally the 
T 3 dependence of 7^ is consistent with the prediction of 
a theoretical analysis |15| of quasiparticle scattering and 
recombination in the presence of antiferromagnetic cor- 
relations. Although the rates were calculated within the 
random phase approximation, the asymptotic T 3 behav- 
ior is a universal feature of quasiparticle scattering in a 
d-wave superconductor |lq| . 

The simplest interpretation of the connection between 
7(0) and quasiparticle dynamics is that AR(t)/R reflects 
the decay of the nonequilibrium quasiparticle density. 
According to this view, 7(0) is the rate at which quasi- 
particle pairs scatter into the condensate and jth is its 
value in thermal equilibrium. Although such an inter- 
pretation appears to be consistent with the prediction of 
Ref. 16, this description of the decay of AR/R is not 



internally consistent. The inconsistency arises because 
quasiparticle interactions alone cannot reduce either the 
excess particle or energy density. For quasiparticles to 
scatter into the condensate they must shed their excess 
energy, which for antiferromagnetic interactions appears 
as a spin fluctuation. However, in the low-energy regime 
(lu << J) the spin fluctuation is essentially the same 
entity as a quasiparticle-hole pair. Therefore such a 're- 
combination' process actually scatters one pair into an- 
other, with no net change in the number or energy den- 
sity. This is analogous to a problem in the interpretation 
of transport experiments, where quasiparticle scattering 
alone cannot reduce the total momentum |lq| . 

The 'true' recombination rate, in the sense of recovery 
of the of equilibrium quasiparticle density, would likely 
differ markedly from the reflectivity decay rates shown in 
Fig. 4. Relaxation of the quasiparticle density requires 
that excess energy flow irreversibly from the electrons to 
the lattice JL7J. We expect that process to be extremely 
slow on the time scale of quasiparticle scattering and to 
exhibit single rather than two-particle kinetics. This ex- 
pectation is based on the physics of the 'phonon bottle- 
neck' in s-wave superconductors 113] , an effect which may 
be even more restrictive in d-wave superconductors. 

The bottleneck arises because the process in which a 
quasiparticle pair emits a phonon and scatters into the 
condensate is reversible — the emitted phonon can im- 
mediately regenerate the quasiparticle pair. The quasi- 
particle density relaxes at the rate at which the typical 
phonon energy decreases below the threshold for pair cre- 
ation. Not only is this rate slower than the quasiparticle 
scattering rate, it is also approximately independent of 
the excess quasiparticle density. 

The relaxation of the quasiparticle density may be even 
slower in the cuprates than in conventional superconduc- 
tors. In s-wave superconductors the phonons decouple 
from the electrons when their energy becomes less than 
2A. In the d-wave superconductors, where there is no 
threshold for pair creation, the typical phonon energy 
must degrade to ~ ksT in order for the electronic system 
to reach equilibrium. The extremely weak coupling be- 
tween nodal quasiparticles and phonons compounds this 
effect. The decay of a quasiparticle pair into an acoustic 
phonon is kinematically forbidden because the quasipar- 
ticle velocity is larger than the sound velocity [17j . 

If recombination is extremely slow, why does AR/R 
decay on a picosecond time scale? The remaining pos- 
sibility is that the decay corresponds to thcrmalization, 
rather than recombination. If AR/ R is not a direct mea- 
sure of the quasiparticle density, but proportional instead 
to the £( 2 ) as we have argued previously, it can decrease 
during thermalization, despite a constant, or even in- 
creasing, density of quasiparticles. Microwave measure- 
ments show that S*- 2 ) for nodal particles is small, because 
their spectral weight is confined to a very narrow (~10 
/ieV) Drude peak Ea]. The observed ~100 meV spectral 



weight shifts 0, and consequently large values of T,^ 2 \ 
are likely due to the antinodal quasiparticles. These con- 
siderations suggest that the decay of AR/R measured 
at 1.5 eV tracks the conversion of quasiparticles from 
a nonthermal distribution, which has substantial T,( 2 ', 
to a quasithermal distribution in which £( 2 ) is a factor 
10 6 smaller. Measurements in which the nonequilibrium 
state is probed in the frequency range of the Drude con- 
ductivity |17| show much longer decay times (typically 
milliseconds), consistent with an experiment which se- 
lectively probes the nodal quasiparticle density. 

In summary, we present measurements of the photoin- 
duced change in reflectivity in an underdoped crystal of 
YBCO, and discuss a new perspective on this effect. We 
interpret the decay rate of AR/R as the inelastic scatter- 
ing rate of quasiparticle pairs, which sets the time scale 
for thermalization of the nonequilibrium population. The 
thermalization rate decreases with the development of 
the pseudogap, becomes excitation density dependent in 
the superconducting state, and varies as T 3 in the low 
excitation limit. Coordinated measurements of photoin- 
duced changes in optical response over a broad spectral 
range, from terahertz to visible frequencies, are needed 
to test this perspective. 
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